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Abstract 

In a recent paper, Filatrella et al. [Phys. Rev. B 75, 54510 (2007)] report results of numerical 
calculations of energy barriers for flux quanta propagation in two-dimensional arrays of Josephson 
junctions with finite self and mutual inductances. To avoid complex numerical calculations, they 
use an approximated inductance model to address the effects of the mutual couplings. Using a full 
inductance matrix model, we show that this approximated model cannot be used to calculate the 
energy barrier of interacting vortices. The authors' method of comparison between numerical data 
and theoretical predictions is also discussed. 
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The static properties of single-vortex states in two-dimensional arrays of Josephson junc- 
tions have been analyzed by Phillips et al. [l] when the self field effects are important. Using 
numerical calculations, they found that the full inductance matrix is necessary to recall the 
magnetic properties of the vortex that are similar to the magnetic properties of a vortex in a 
thin superconducting film. They also showed that it is sufficient to include self and nearest- 
neighbour terms in the inductance matrix to model the energy barrier for the cell-to-cell 
motion of an isolated vortex. In fact depending on the structure of the vortex near the core, 
the energy barrier is less affected by the mutual inductances, while the vortex properties 
that are related to longer length scale, such as the vortex shape or the current distribution, 
are more strongly affected. Moreover, they stressed that when a vortex model includes only 
the self terms of the inductance matrix, the magnetic properties of a vortex are similar to 
those of a vortex in a bulk superconductor. This vortex model is sometimes referred to as 
the Nakajima-Sawada model 2|. 

In a recent paper Filatrella et al. 3|] present a numerical study on the effect of noise on 
vortex propagation in two-dimensional arrays of Josephson junctions. They evaluate the 
energy barrier in the vortex-(anti)vortex interaction. Two models of vortex are considered: 
the Nakajima-Sawada (NS) model and the mutual coupling (MC) model, which is an ap- 
proximated inductance model, developed by the authors, to address the effects of the mutual 
couplings between the array elements and thus avoid the complex numerical calculations of 
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a full inductance matrix model [4]. They calculate the energy barriers for the cell-to-cell 
motion of an isolated vortex to test the goodness of the MC model. The compatibility of 
the findings with the numerical results of Phillips et al.pj] is considered by the authors as a 
good test for the MC model. 

We observe that although a simple self and nearest-neighbour inductance model is suffi- 
cient for consistency with the results of Phillips et al.|l] it is necessary to demonstrate that 
an approximate inductance model, such as the MC model, is equivalent to a fully inductance 
model [sl. 

Figure 5(d) of Filatrella et al. shows the energy barrier for the cell-to-cell vortex motion 
as a function of the distance between two vortices. This is the case in which the energy barrier 
adds the contribution of the vortex-vortex interaction. The comparison between the data 
obtained from the NS model with those from the MC model predicts that the energy barrier 
of two interacting vortices in a bulk superconductor cannot be distinguished from that of 
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m. In other words, the characteristic 



does not produce a clear signal for 



two interacting vortices in a thin superconducting fi 
long-range interaction of two vortices in a thin film 
the vortex motion. 

To show that this result depends on the method developed by the authors to address 
the mutual magnetic couplings we have repeated the analysis using a full inductance 

matrix model [7| . 

In Figured] we show the comparison between the NS model with the full inductance model 
at (3l = 0.5. The data clearly show that the energy barrier of two interacting vortices in 
a full inductance matrix model is different from that of two interacting vortices in a model 
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that includes only the self terms of the inductance matrix. Unlike Filatrella et al.[3|, we 
find that the minimum distance for a static pair of vortices also depends on the model, as 
consequence of the different behaviour of the interaction. This distance is equal to three 
cells for the NS model and to four cells for the full inductance matrix model. 

Moreover, from another point of view, the force of interaction between two vortices de- 
)ends on the vortex current produced by the first vortex at the centre of the second vortex 
Sj and, as Phillips et al.[li| have shown, the correct calculation of the current distribution 
requires to include all terms in the mutual inductance matrix. Therefore, an approximated 
inductance model is not able to simulate the long-range interaction between vortices. 

The second consideration on the study presented by Filatrella et al. ^ is the comparison 
between numerical data and theoretical predictions. 

For a single vortex-state the energy barrier height Aeo is given by the energy difference 
between a vortex centered on a i unction and a vortex centered in the middle of the cell to 
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which the junction belongs [9|. 

For two interacting vortices at a relative cell distance d, the energy barrier height is 
calculated, by the authors, as the energy necessary to shift one of the vortices by half a cell 
^ and it is splitted in two terms: 

Ae = Aeo + Ae^"* (1) 

where Aeo is the term due to the vortex interaction with the network and Ae*"* the term 
due to the interaction between vortices, connected to the vortex- vortex interaction potential 
e*"*. In particular, regarding the MC model, Filatrella et al.fsj analyze the behaviour of 
Ae*"* by fitting the numerical data with a power law as a function of the distance d between 
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the vortices: 

- Ae'"* = 4 (2) 
and find A ~ 3, in units of Josephson energy, and a ~ 3.1. Tlie last value is compared witli 



10| for the vortex- vortex 



the theoretical prediction a*^ = 1 obtained by Dominguez and Jose 
interaction energy of a full inductance model. 

We observe that this comparison is not correct because the two quantities are not phy- 
sically equivalent, since the prediction by Dominguez and Jose jl^ refers to the interaction 
energy, i.e. the energy necessary to put two vortices, starting from infinity, at a finite 
distance d: 

Ae^"* = e*"*(rf) - e™*(oo) = ^ (3) 



while Ae**^* for Filatrella et al.[3| is the contribution of the vortex- vortex interaction to the 
energy barrier, i.e. the change of the interaction potential when, starting with two vortices 
at a finite distance rf, one of them is displaced, toward (d ^ d — 1/2) or apart [d (i+ 1/2), 
by half a cell through the lattice barrier. Therefore, the right predictions for the two cases 
are: 

- A6^"* = e^"*(rf) - e^"*(^ - 1/2) = ^ - , ^ , d^d-l/2 (4) 

d d 1/2 

-Ae*"* = e*"*(rf) - e'^'Hd + 1/2) = 4 - , ^ , d^d+ 1/2 (5) 

d ct + 1/2 

(6) 

Using a full inductance matrix model 7| we have analyzed the behaviour of Ae*"* when one 
vortex is moved apart by half a cell in a 75 x 75 square array at I3l = 0.5. Fitting the data 
reported in Figure [2] with the function: 

A A 

_ ^ int ^ /7^ 

rf" d" + l/2 ^' 

we find A=7.87 units of Josephson energy and a = 0.85. The result is in good agreement 
with the theoretical prediction by Dominguez and Jose |lol |. 

This work has been partially supported by MIUR PRIN 2006 under the project: "Macro- 
scopic Quantum Systems - Fundamental Aspects and Applications of Non-conventional 



4 



Josephson Structures" 



[1] J.R. Phillips, H.S.J, van der Zant, J. White, and T.P. Orlando, Phys. Rev. B 47, 5219 (1993). 
[2] K. Nakajima and Y. Sawada, J. Appl. Phys. 52, 5732 (1981). 
[3] G. Filatrella, S. Girotti, and G. Rotoh, Phys. Rev. B 75, 054510 (2007). 
[4] G. Rotoh, S. Girotti, and G. Filatrella, IEEE Trans. App. Supercond., vol. 17, Issue 2, part 
3 (2007), 353. 

[5] C. De Leo, private communication to G.Filatrella and G. Rotoli, January 2006. 
[6] J. Pearl, Appl. Phys. Lett. 5, 65 (1964). 

[7] C. De Leo, IEEE Trans. App. Supercond., vol. 18, Issue 4, (2008), 1769. 
[8] V. V. Shmidt and G. S. Mkrtchyan, Sov. Phys. Usp. 17,170 (1974). 
[9] C.J.Lobb,D.W.Abraham,and M.Tinkham, Phys. Rev. B 27, 150 (1983). 
[10] D. Dominguez and J.V. Jose, Phys. Rev. B 53, 11692 (1996). 



5 



FIGURE CAPTIONS 



FIG. 1: Energy barriers for two vortices in units of Josephson energy as a function of the distance 
between the vortices. Comparison of the NS model (filled symbols) with the full inductance matrix 
model (empty symbols) for (3l = 0.5 in the case of vortex-vortex (circles) and vortex-antivortex 
(triangles) interactions. 

FIG. 2: Normalized interaction energy contributes to the pinning barriers for two vortices. Vortex- 
vortex pair in a full inductance matrix model with (3l = 0.5. The solid line has been obtained 
fitting the data with the function of Eq. [7| 
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